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Abstract. An analytically calculable model ^ has been proposed recently for 
the simultaneous description of particle spectra and correlations, which is able 
^ . to explain in a natural manner the mj scaling of the radius parameters of the 

Bose-Einstein correlation functions. We present here preliminary results on 
. fitting this model to data sampled by NA44 collaboration in central S + Pb 

' reactions at CERN SPS, which indicate that both pions and kaons freeze-out 

at a surprisingly low temperature of T ~ 108 MeV. 

^ ■ 1. Introduction 

Recently there has been much interest in the measurement and the calculation of 
^ i' the Bose-Einstein correlation functions (BECFs) and those of the invariant momen- 

O . tum distributions (IMDs) for rapidly expanding systems with flow and temperature 

profiles. The expansion may result in strong correlations between space-time and 
momentum space variables implying that the Bose-Einstein correlations are in gen- 
eral not measuring the geometrical sizes of big and expanding finite systems. The 
geometrical sizes play an important role in the description of the invariant momen- 
' tum distributions. That is why a simultaneous fit of IMDs and BECFs gives good 

constraints on the parameters of the model. 

We will discuss an attempt to fit simultaneously the data of S-f Pb interac- 
tions analysed by the NA44 collaboration They showed an unexpected, scaling 
behaviour for the radii parameters of the BECF, Rside: Rout and Riong- These pa- 
rameters turned out to be proportional to 1/ y^mj, where mt is the mean transverse 
mass of the identical particle pair. The same proportionality factor was valid for 
both pions and kaons. The NA44 collaboration has also measured the IMDs for 
the same S-|-Pb interactions using similar trigger conditions. Data are not yet fully 
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normalized but that is not so serious for our purpose as we will only fit the shape 
of the IMD anyway. We have two reasons to only take the shape into account: 

(i) The slightly different trigger conditions for data sampled to measure IMD 
and for BECF might give slightly different normalizations but the shape is observed 
to be uneffected 

{a) The expansion model is applicable only to the central finite core creating 
particles, not to all particles measured. This core is characterized by hydrodynamic, 
collective behaviour and small lengths of homogeneity. This is in contrast to the 
halo, created by the decay products of long-lived resonances, which is characterized 
by a free streaming and decay of these resonances and correspondingly with large 
lengths of homogeneity proportional to the decay time of the long-lived resonances, 
[Q . The A parameter which measures the strength of the correlation gives a handle 
to estimate the fraction of particles which are produced in the core |, |. When A 
is independent of rrit the IMD of the core particles is proportional to the IMD of 
all particles measured. The S+Pb data show an approximately nit independent A 
parameter, A7r+ = 0.56 ±0.02 and 0.55 ±0.02 at the quite different mean transverse 
momenta of 150 MeV and 450 MeV, respectively. 

We will shortly present the model of the core, present the results of a fit to 
three sets of IMDs and the corresponding radii measurements using (i) an analytical 
approximation and (ii) a numerical integration method. We will summarize with 
some conclusions. 



For central heavy ion collisions at high energies the beam or z axis becomes a 
symmetry axis. Since the initial state of the reaction is axially symmetric and 
the equations of motion do not break this pattern, the final state must be axially 
symmetric, too. However, in order to generate the thermal length scales in the 
transverse directions, the flow field must be either three-dimensional, or the tem- 
perature must have significant gradients in the transverse directions. Furthermore, 
the local temperature may change during the the duration of the particle emission 
either because of the re- heating of the system caused by the hadronization and/or 
intensive rescattering processes or the local temperature may decrease because of 
the expansion and the emission of the most energetic particles from the interac- 
tion region. We consider the following emission function for high energy heavy ion 
reactions as a model of the core: 



2. A Model of the Core 




9 



(1) 




where (i^S''(a;)pp describes the flux of particles through a finite layer of freeze-out 
hypersurfaces. We assume that this layer is parameterized by r = ^/t"^ — where 
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the random variable t is characterized by a probabihty distribution, such that 

(i'*S''(a;)p^ — nit cosh[rj ~ y] H{T)dT tqcIt] dr^ dvy , (2) 

where the finite duration is modelled by H{t) cx exp(— (r — to)^/(2At^)). Here 
To is the mean emission time, At is the duration of the emission in (proper) time. 
The four-velocity and the local temperature and density profile of the expanding 
matter is given by 

u{x) — I cosh[77] cosh[77t], sinh[77t] — , sinh[77t] — , sinh[?7] cosh[77t] I , (3) 
\ n rt J 




T{x) To \ 2t^ J \ j ' To 

T{x) Tc 



(5) 



where iJ.{x) is the chemical potential and T{x) is the local temperature character- 
izing the particle emission. 

Note that in Ref. [Q] we have considered the Boltzmann approximation to 
the emission function given by the present eqs. (|l]- |^). Thus the results obtained 
in Ref. |^ correspond to Boltzmann approximations. Here we indicate some new 
analytic results for the single particle inclusive spectra which are obtained without 
the full Boltzmann approximation. 

The emission function of eqs. (0 - ^ corresponds to a longitudinally expanding, 
finite system, with transverse and temporal temperature profile and with a linear 
transverse flow velocity profile. The decrease of the temperature in the transverse 
direction is controlled by the parameter a, while the strength of the transverse 
flow is controlled by the parameter b. The parameter d controls the strength of 
the change of the local temperature during the course of particle emission. This 
parameterization corresponds to a Gaussian approximation of the local flow, inverse 
temperature and chemical potential profiles, as discussed in detail in Ref. [Q. For 
the case of a = d = we recover a model also discussed in Ref. ^] , while for the 
case of d = the model of Ref. 0| is reproduced, while for & = d = we obtain a 
version of a model proposed in Ref. Q . 

Using a saddle point approximation similar to the one presented in Ref. 
but keeping the Bose-Einstein form of the emission function, we find that the IMD 
of the particles produced from the core plus the halo can be expressed as: 

d^n 1 d^Uc 



dy drn^ JY^ dy dm^ ' 



(6) 

^ '^'"^ - 5 , ^ . (7) 



TTdydrUt (2n) / fi{xs) 

'''''' i Tixs) Tixs) 
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This result can be considered new as compared to the ones presented in Ref. 0], 
but applying a Boltzmann approximation to this result reproduces the expressions 
given in Ref [Q. Parts of this expression are explained piece by piece below. The 
intercept parameter A», which relates the total number of particles to the number 
of particles originating from the core, is the effective intercept parameter of the 
Bose-Einstein correlation function. The intercept parameter A* is not prescribed 
by the model, it has to be taken from the correlation measurements. Note that A* 
is in general rapidity and transverse mass dependent and thus may depend on the 
kinematic region where the data are determined. The effective volume factor V^, [xs) 
measures the invariant volume element of the particle emission at the time of the 
last interactions. This factor can he very strongly mt dependent, and its possible 
mt dependence should he taken into account in the future data analysis M. This 
invariant, effective volume factor is related to the HBT radius parameters |1[ as 

V;(x.) = (27rro2Ary^)^ (2^i?2) At._ 

Ar 

The effective source dimensions in the transverse, longitudinal and temporal direc- 
tions were obtained in Boltzmann approximation in Ref. [nj as 



1 1 1 



It^ JXq it' 

1 1 



—2" cosh[77s], (9) 



Ayy^ Arf At]'^ cosh[77, 

1 1 1 



cosh[77,] -—2, (10) 

cosh[r/3] , (11) 



- 2 y i.2~' (12) 



Ar^ Ar2 Ar| 

a^ + b'' mt 



Ar?J = ^, and Ar^ = 4 — - (13) 
mt d mt 



while the Cooper-Frye flux factor contains n'^, the normal-pointing unit vector of 
the hypersurface family c?''S'^(x). This Cooper-Frye flux factor is given by 

p'^n^{xs) = mfCOsh[77s] (14) 

In the present status report we do not yet consider the effects of Bose-Einstein 
statistics on the saddle-point and on the radius parameters of the Bose-Einstein 
correlation functions. In Boltzmann approximation, the saddle point x^ is charac- 
terized by Ts = To, ris = (j/o - J/)/[l + A.jf{mt/To - 1)], r^.s = PtbRl{r]s)mt / {tqTq) 
and Tj, s = 0, if the saddle-point equations can be linearized in the LCMS |^ around 
r^ — 0. For the NA44 experiment, the mid-rapidity is located at j/o = 3. 



Hydrodynamic Parameterization of NA44 Data 



225 



The same length scales, which determine the effective, invariant volume factor, 
K, shall determine the radius parameters of the BECF. These length scales, indexed 
by subscript *, are the homogeneity lengths in the considered spatial or temporal 
directions. From the above expressions it is clear that they are determined by 
an interplay of the geometrical scales, Rg, A77, and Ar and the so-called thermal 
scales, i?T, A77T, and Att, which are indexed with subscript T and stem from the 
Boltzmann factor eyi]i[p^u^{xs) lT{xs)\- In Boltzmann approximation, utilizing a 
solution of the linearized saddle point equations one obtains in the mid-rapidity 
region, for y = yg, in the LCMS system 

Ride = Rl (15) 

Rlut - i?^+/3?Ar2, (16) 
Rl - T-2Ar,2 . (17) 

The BECF can be expressed as 

C{Q- K)^l + K eM-RldeQlde - RlutQlut - RlQl - '^Rlm.LQoutQL) (18) 

For the NA44 data, which are sampled in the central region, y w 3, we have 
R^ut i ~ 0. A more detailed discussion on the HBT radii was given in Ref. ||l|. 

From the above expressions one can figure out that the radius parameters in 
general may have complicated dependence on the rapidity and on the transverse 
mass of the particle pairs. This is partly due to the rapidity and transverse mass de- 
pendence of the = ^LCMS variable, the cross-term generating hyperbolic mixing 
angle ||l| defined in the LCMS frame, and partly due to the simultaneous presence 
of mt dependent and mt independent terms in eqs. (|l2| - |l^). We see also, that 
all radii parameters may be proportional to l/y/rrit in the case when the total size 
and time spread of the particle emission are much larger than the corresponding 
thermal sizes, i?T, At^t and Att, respectively. 

There is no a priori reason to assume that STr{x,p) = Sk{x,p) i.e. that the 
emission functions have the same form for pions (tt) and kaons (K). Freeze-out 
arguments are based on the estimated value of the hadronic cross-sections and the 
pions have larger cross-sections than the kaons. On this basis, one would expect 
that pions will freeze-out later, from a larger volume, than the kaons. 

Our results on the parameters of the two-particle correlation function indi- 
cate, that even if the geometrical radii for the pion emission function and the kaon 
emission function were different, they may play only a minor role if they are suf- 
ficiently large as compared to the thermal length-scales. Thus the parameters of 
the Bose-Einstein correlation function may obey an rrit scaling even if the emission 
funtions of the various hadronic pieces were different [|l^ . 

The radii parameters are mainly sensitive to the smallest of the 2 scales 
present, in contrast to the IMDs which are sensitive to the biggest of the scales. The 
model is discussed in detail in Ref. [Q, where the spectrum has been calculated 
in a Boltzmann approximation, allowing for a detailed analysis of the influence 
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of the thermal and geometrical length scales on the rapidity width and the slope 
parameter of the momentum distribution. 

3. Fitting Data with Analytical Formulas and Using Numerical 
Integrations 

In the previous section, we have discussed some new features of the Bose-Einstein 
statistics on the momentum distribution obtained in the saddle-point approxima- 
tion. In this section, we present a status report about fitting the Boltzmann ap- 
proximations to eq. (|l|) to NA44 data, both analytically and numerically. Our 
fitted formulas thus are not taken from the previous chapter, but from Ref. 
which corresponds to an improved saddle-point approximation to the Boltzmann- 
approximated emission function. Work is in progress to study the effects of Bose- 
Einstein statistics on the analytic and numerical fitting of these NA44 data. In a 
qualitatively similar analysis of the NA49 data the effects arizing from the Bose- 
Einstein statistics were reported to be significantly influencing the quality of the fit 

Thus we have tested the model of ref. [y by a x^-?it to the NA44 data on 
the parameters of the kaon and pion BECF as well as on the kaon and pion IMD. 
In the fit, the data points for the kaon and pion spectrum did not include yet the 
systematic errors, which were first estimated in Refs. We have fitted the 

analytical results for the spectrum obtained in a Boltzmann approximation using 
the modified saddle point method, as described in Ref. Q. 

Our results are still preliminary because of the following reasons: 

i) PiX the time of the fitting the systematic errors on the pion and kaon IMD were 

not yet available, which resulted in underestimated errors on the spectrum. 
The small errors in turn imposed a too strong requirement on the fitting, and 
increased the relative weight of the spectrum data points with respect to the 
radius parameters, for which both the systematic and the statistical errors 
were determined experimentally. As a consequence, the preliminary results 
discussed below underestimate the errors of the fitted parameters. 

ii) We are currently developing a method which supposedly will allow us to analyt- 

ically estimate the errors made when the saddle-point integrations are carried 
out. This could be potentially very helpful to restrict the domain of the pa- 
rameters in the analytical fitting procedure to that part of the parameter 
space where the analytical results correspond to the numerically integrated 
distributions |]l3| . 

in) We have note yet included the effects of the Bose-Einstein statistics into the 
fitting, and finally 

iv) We currently realized that the hypersurface has a rather complicated shape 
when comparing our models to NA44 pion and kaon data. We hope that 
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the distinction among the various almost degenerate but physically different 
minima can be improved in the near future by including data ondata on 
heavier particles into the data analysis. 

The final results of the data analysis, to be obtained by eliminating the problems 
in items i) - iv), will be reported elsewhere. 

With the reservations as listed above, we have found the parameter values, 
summarized in Table 1, by applying the CERN MINUIT multidimensional fitting 
algorithm to determine the lowest values of /NDF in the parameter space of the 
model. When we fit to pion data separately or only to kaon data, separately, wc do 
not find a fit which were statistically better than the simultaneous fit to both the 
pion and kaon sample, see Table 2. The overall 'x^ /NDF does not change. 

Table 1. Fitting HBT radius parameters and spectra from the model with 
analytical approximations, where the domain of the validity of the approxi- 
mations is not taken into account, and with the numerically integrated radius 
parameters and spectra. 





Analytical 
(preliminary errors) 


Numerical 
(preliminary errors) 


x'/ndf, full fit 

Ax^ /bins, low pt 
Ax^ /bins, high pt 
Ax^ /bins, kaon 
Ay/^ /bins, radii 


207/126 

77/42 
78/35 
41/51 
12/9 


237/126 

104/42 
76/35 
40/51 
18/9 


To [MeV] 


140 ± 2 


108 ± 10 


TO [fm/c] 


4.6 ± 0.2 


7.0 ± 0.5 


Rg [fm] 


7.0 ± 1.1 


4.4 ± 0.2 


At] 


2.5 ± 1.3 


1.5 ± 0.6 


At [fm/c] 


3.2± 1 


2.4 ± 1.3 


/ 


0.89 ± 0.03 


0.88 ±0.12 




0.52 ± 0.08 


0.68 ±0.14 


rf2 


499 ± 50 


130 ± 8 



In Fig. 1 we show the result of the fit with the analytical expressions, displayed 
by a dashed line. The measured data points together with their error bars are 
also shown. Parameters a and b appear in a combined manner in the analytic 
expressions for the invariant momentum distribution, that is why the new variable 
/ = P/ia"^ + 6^) is introduced in the fit. 

The analytical formulas give a slightly better adjustment to the data than the 
description obtained by the numerically integrated emission function. This may 
either indicate that the analytical fit develops a fake minimum, or alternatively 
that some physics appears effectively in the analytic approximations which was 
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still missing from the emission function of eq. (|l|). The analytical formulas thus 
correspond to the source given by eq. (|^) under certain restrictions only. A detailed 
numerical study indicated that the conditions 

= ^5-7 — 9- < 0.6 

To T^A-q^ 

and Ari^{y,mt) < 0.9 have to be satisfied simultaneously in order to reach less 
than 20 % relative error in the HBT radius parameters and the slope parameter of 
the IMD. This condition is violated by our analytic fit to the pion sample: from 
< 0.6 we obtain mt,7r < 200 MeV, while from A?7* < 0.9 one obtains nit^Tr > 300 
MeV, and these conditions cannot be satisfied simultaneously. The same conditions 
imply mt.K < 600 MeV, which is satisfied for the NA44 kaon data. 

Table 2. Fitting either NA44 n spectrum and radius parameters or K radius 
parameters and spectrum, with the numerically integrated model, in Boltz- 
mann approximation. 



Only TT Only K 

(preliminary errors) (preliminary errors) 



xVnd/, full fit 

/bins, low pt 
Ax^/bins, high pt 
Ax^ /bins, kaon 
Ax^ /bins, radii 


195/76 


40/48 


103/42 




75/35 






40/51 


16/6 


0/3 


To [MeV] 


112 ± 2 


133 ± 6 


To [fm/c] 


6.6 ± 0.3 


6.4 ± 0.4 


Rg [fm] 
Arf 


4.1 ± 0.1 


12 ± 1 


=2.2 


=2.2 


At2 [fm2] 


=5.5 


=5.5 


/ 


1.0 ±0.2 


0.52 ± 0.1 




0.532 ±0.02 


0.58 ± 0.02 




=130 


=130 



Note that the restrictions given by the IMD data points, which do not yet 
have the systematic errors, seem to be so strong that the HBT mj-dependent radius 
parameters from the fit are forced to deviate from the Rl — Rside — Rout oc l/^/roJ 
rule, which is well satisfied by the NA44 data points. If we do tried not to fit the 
IMD of pions and kaons, our model were able to perfectly fit the 9 NA44 HBT data 
points, with large allowed region for the geometrical source sizes. 

On this basis we expect that the inclusion of the systematic errors for the IMD 
of pions and kaons will improve the agreement between the model and the data. 
Note, however, that a fit to the measured particle spectrum is rarely obtained in 
such a good quality as obtained here. The fit with numerical integration is shown 
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on Fig. 1 with a full line. The of the fit is somewhat high but still acceptable. 
The fit with the analytical formulas is shown on Fig. 1 with a dashed line. The 
difference between the two fits is only seen in the low region for the pions. 




Pj /GeV/ 

Fig. 1. The upper half shows the two pion IMD, the lower half the kaon 
IMD. Data are here displayed versus pt and they are unnormalized. The full 
line shows the result of the fit with numerical integration, the dashed line 
the result using the analytical formulas. The rapidity of the data points is 
not shown. The acceptance is restricted to a small region in rapidity which 
changes with pt specially for the pion data More than one data point for 
a Pi means neighbouring rapidity bins. The rapidity of the pions varies from 
4.0 to 2.6, the rapidity of kaons spans 3.4 to 3.1. Data points at the edge of 
acceptance are not used. 
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0.6 U.» 1 

m,-in, /GeV/ 



0.6 0.8 1 

m,-m, /GeV/ 



Fig. 2. Full line shows the fit using the numerical integration, dashed line 
shows the result of the analytical formulas with the parameters of the same fit. 
Data points are the radii measurements of NA44 |^. Unfilled circles are used 
for pion data, filled for kaon data. 



Fig. 2 shows the fit with numerical integration to the radii parameters and 
also the difference between the numerical integration and the analytical calculation 
using the parameters of the fit with numerical calculation. We see important differ- 
ences mainly for the IMD. The radius parameters are reproduced by the analytic 
expressions with less than 20 % relative error at the transverse mass of the data 
points, and the slope parameters of the IMD for the pions is also reproduced up to 
a 20 % relative error. 
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When we fit the analytical formulas without restricting their domain of valid- 
ity, the results are given in Table 1. Unfortunately, those parameter values which 
correspond to the best description of the data by the analytical formulas are outside 
of that region of the parameter space, where the analytic approximations correspond 
to the numerically evaluated emission function. Either the model emission function 
needs to be further adjusted in such a manner that its numerically integrated spec- 
trum and HBT radius parameters will closely correspond to the presently found 
analytical formulas, or, if the emission function is fixed at its present form, one 
should not allow for fitting the data with the analytical formulas in an unrestricted 
domain of the parameter space. A more detailed analysis is needed, which can take 
into account that the errors of the analytic formulas increase strongly if we leave 
the domain of the 20% relative errors. This domain is given by the requirements 
^ < 0.6 and Ar/^ < 0.9, see Ref. @ for further details. 

The fit implies that the model parameters a, b and d take values significantly 
different from 0, approximately 0.3, 0.8 and 11, respectively. 

4. Conclusions 

This hydrodynamical model gives a fair description of data. 

The analytical approximations provide less than 20% relative error for the 
HBT radius parameters and the slope parameters of the spectra in that part of the 
parameter space where the NA44 data are. 

If we fit the approximative analytic results to data without estimating the 
errors of the approximations in the region, where the fit is performed, we find a 
fake minimum, where the extrapolated analytic results give a very good description 
of the data. However, at this part of the parameter space, no minimum is found 
when numerical integration is invoked. We find no evidence for different model 
parameters for pions or for kaons. 

In a preliminary analysis, we find a surprisingly low central temperature at 
the mean time of the last interactions, Tq = 108 ± 10 MeV. This is in agreement 
with our preliminary results presented at Quark Matter'95 0. Note that a similar 
low value of Tg = 92.9 ± 4.4 MeV has been found at AGS energies from a detailed 
analysis of the BECF and the spectrum of both pions and kaons |Q. However, 
we would like to be careful and not to draw very strong conclusions untill we can 
find a very well defined minimum by adding more data, utilizing the experimental 
estimate of the systematic errors on the spectrum, and untill we have included the 
effects of quantum statistics on the single-particle spectra into the data analysis. 
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